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Abstract 

We study electronic structures of two-dimensional quantum dots in high magnetic fields using the density-functional theory 
(DFT) and the exact diagonalization (ED). With increasing magnetic field, beyond the formation of the totally spin-polarized 
maximum density droplet (MDD) state, the DFT gives the ground-state total angular momentum as a continuous function 
with well-defined plateaus. The plateaus agree well with the magic angular momenta of the ED calculation. By constructing 
a conditional wave function from the Kohn-Sham states we show that vortices enter the quantum dot one-by-one at the 
transition to the state with the adjacent magic angular momentum. Vortices are also observed outside the high-density 
region of the quantum dot. These findings are compared to the ED results and we report a significant agreement. We study 
also interpretations and limitations of the density functional approach in these calculations. 
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Quasi-two-dimensional (2D) quantum dots fabri- 
cated into semiconductor heterostructures are ex- 
pected to be basic components in the future nano- 
electronics [1]. Moreover, they provide a unique play- 
ground in physics with a rich diversity of quantum 
effects among interacting electrons [2]. The reason is 
especially the strong influence of magnetic fields on 
the electron states. Already for fields attainable in the 
laboratory environment the interaction with electrons 
is of the same order as the electron-electron inter- 
actions. The understanding of experimental findings 
in quantum dots and the prediction of new phenom- 
ena have inspired a huge amount of theoretical and 
computational work. 

Recent developments in the electronic structure the- 
ory of 2D quantum dots in magnetic fields have shown 
a rich variety of phenomena related to vortices [3,4,5,6] . 
These vortices correspond to magnetic field quanta 



with rotating currents of charge around them and a 
zero in the wave function at the vortex centra. By go- 
ing around one vortex the wave function gains a phase 
change of 27m, where n is the winding number. In the 
fractional quantum Hall effect (FQHE) the state of the 
matter is described by Laughlin wave functions which 
attach vortex zeros at each electron. Electronic struc- 
ture calculations of charge droplets in quantum dots 
have predicted the appearance of additional free vor- 
tices which cluster near electrons [3,4]. 

In this paper we perform a detailed study of vortex 
formation in a circulary symmetric 2D quantum dots 
in a parabolic external confinement V c (r) — |ra*a;or 2 . 
The model the system by an effective- mass Hamilto- 
nian 
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where N is the number of electrons in the dot, A is 
the vector potential of the external perpendicular mag- 
netic field B, m* the effective electron mass, and e is 
the dielectric constant of the medium. We use the ma- 
terial parameters of GaAs, m*/m e = 0.067 and e/eo = 
12.4. In high magnetic fields the quantum dot is spin- 
polarized due to the Zeeman effect and the maximum 
density droplet state (MDD) is formed [7]. This state is 
a finite- size precursor of the v = 1 quantum Hall state. 
In parabolic external confining potential the state is 
formed by the lowest- Landau- level (LLL) orbitals with 
angular momenta 1 = 0,— 1,...,— iV + 1, and the total 
angular momentum L equals to —N(N — l)/2. For a 
stronger B, ED has shown ground states to occur only 
at certain "magic" L values, and L shows a stepwise 
structure as a function of B [8,9,10]. 

We solve for the ground state electronic structure 
using the spin-density-functional theory (SDFT), the 
current-density- functional theory (CDFT) [11], and 
the exact diagonalization method (ED). The density- 
functional approaches are implemented in the real- 
space without symmetry restrictions [12]. The main 
results were calculated using 128 x 128 grid points 
but similar results were obtained using larger grids up 
to 256 x 256 grid points. High spatial resolution of 
fine grids is necessary to accurately describe systems 
with several vortices inside the quantum dot. The 
exchange- correlation effects are taken into account 
using local approximations for the spin densities and 
the vorticity in the 2D electron gas [13,14]. 

The Hamiltonian of the system is rotationally sym- 
metric and therefore the particle density, defined as 
|^(r)| 2 , should also be rotationally symmetric. Calcu- 
lations, however, show that the particle density in the 
density-functional formalism is not necessarily sym- 
metric. In this work we adopt the approach that even 
if broken symmetry solutions cannot be identified as 
the true particle density we may use them to find out 
correlations in the system. Otherwise, the correlations 
may remain hidden in the one-particle picture of the 
density functional approach. On the other hand one, 
should be cautious with the interpretations since it has 
been shown that the symmetry breaking solutions may 
be wrong mixtures of the eigenstates of the system and 
therefore the interpretation of the solution may be un- 
physical [15]. We therefore emphasize that solutions of 
the DFT should be used as guidelines and approxima- 
tions and the final proofs should be done using exact 




1 ° 10 15 20 

B[T] 

Fig. 1. Angular momentum L of the six-electron QD from 
SDFT (open circles). The plateaus are characterized by 
vortex-holes in the electron density. The number next to a 
plateau gives the number of vortices inside the QD. The 
horizontal lines correspond to ground state magic L values 
from ED. Solutions between the plateaus of nearly constant 
L are charge-density- wave like solutions with fractional L. 

many-body techniques. 

We examine the 6-electron quantum dot (N = 6). 
For an easy comparison with the low magnetic field 
results in Ref. [12] we set the confinement strength to 
coo = 5 meV. The maximum density droplet (MDD) 
state is characterized by a relatively constant electron 
density and currents concentrated at the edge of the 
dot. The transition to the MDD state occurs at 4.8 T. 

The total angular momentum in the SDFT in the 
beyond-MDD domain is shown in Fig. 1. The curve has 
plateaus of nearly constant L. The plateau regions are 
characterized by electron densities with vortex depres- 
sions inside the quantum dot and outside the dot in 
the low density region. The number of vortices inside 
the quantum dot increases by one between adjacent 
plateaus. The plateaus correspond well to the magic 
angular momenta from the ED calculations. In the ED 
four- vortex configuration has the pentagon or hexagon 
symmetry which correspond to angular momenta L — 
—35 and L = —39, respectively. In the DFT solutions, 
however, only the L = — 35 plateau is visible. 

In the density-functional solutions vortices appear 
as static holes which break the rotational symmetry 
of the particle density. The wave function must vanish 
at the vortex position in order to keep the kinetic en- 
ergy finite. In the beyond MDD-domain there are also 
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charge- density- wave (CDW)-like solutions with frac- 
tional L values in the transition regions between the 
adjacent magic angular momentum values. For exam- 
ple, solutions from B ~ 8.3 T to B ~ 11 T have 6 
density maxima at the edge of the dot and a minimum 
at the dot center (see Fig. 2) 

We do a detailed analysis of the vortex solutions in 
the SDFT by constructing a conditional wave function 
of the Kohn-Sham states. We do it iteratively by fixing 
N — 1 electrons to positions {r*}^L 2 an d then calculate 
a determinant of the N occupied Kohn-Sham states 
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The electron positions are moved until the total prob- 
ability of the <3>ks obtains a maximum value. Strictly 
speaking <E>Ks(r) constructed in this way is only an 
auxiliar function which emulates the true conditional 
many-body wave function [3] which is unknown in the 
density functional approach. 

^Ks(r) for the MDD-state correctly assigns one vor- 
tex for each electron. The evolution of this state as the 
magnetic field is increased is shown in Fig. 3. It de- 
picts a vortex hole entering the quantum dot as the 
magnetic field is increased. The CDW solutions can be 
analyzed further by projecting the Kohn-Sham states 
on the Fock-Darwin states [2] and comparing the re- 
sults with those from the exact diagonalization (Fig. 
2). From this data we can conclude that the charge- 
density-wave originates from the combinations of the 
L — — 15 and L — —21 states. The mixture of these 
two states result in a density with six peaks in a form 
of a hexagon. The linear behaviours between plateaus 
can be thought to be unphysical since the exact ground 
state jumps directly from the MDD to L — — 21. How- 
ever the DFT gives the same qualitative picture than 
the ED where the (excited state) solutions between 
L — — 15 and —21 show a vortex-hole moving from the 
outer edge toward the center [16]. The other CDW- 
regions show similar solutions, transport of one addi- 
tional vortex-hole to the center of the dot. 

The triple- vortex solution at B — 15.5 T is shown 
in Fig. 4. The DFT localizes the vortices in the elec- 
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Fig. 2. Charge-density-wave (CDW) solution in a 
six-electron quantum dot at 9.7 T. The total angular mo- 
mentum L « —18.230. Left panel: the electron density. 
Right panel: the projections of the angular momentum to 
Fock-Darwin states for the CDW solution (black bars) are 
compared to the L = — 18 solution from the exact diago- 
nalization method (grey bars). 
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Fig. 3. Conditional wave function from the Kohn-Sham 
states <^Ks( r ) f° r the MDD-state, two charge-density- wave 
states, and a single- vortex state. The contour plots are loga- 
rithmic conditional electron densities, the shadowings show 
the phases of the wave-function, and the crosses mark the 
most probable electron positions. At the lines where shad- 
owing changes from the darkest grey to the lightest grey 
the phase changes by 2tt. 

tron density. The conditional wave function shows vor- 
tices at exactly these positions. External vortices ap- 
pear both in the conditional wave function and in the 
logarithmic plot of the electron density (not shown). 
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Fig. 4. Electron density (grey scale) and the physical cur- 
rent density (vectors) for the CDFT triple-vortex solu- 
tion at 15.5 T (left). Conditional wave function from the 
Kohn-Sham states <E>Ks( r ) f° r the triple- vortex solution 
(right). The pictures have not the same scale. 




Fig. 5. Conditional ED wave function of L = —45 state 
(left) and the $ K s(r) for the L ~ -45.129 state at 21 T 
(right). 

The comparison with the ED results for L — — 30 in 
Ref. [3] shows that the conditional DFT wave func- 
tion given by Eq. (2) correctly predicts the existence 
of three vortices of winding number one inside the ring 
of the fixed electrons. In the ED the vortices start to 
cluster around the fixed electrons for larger angular 
momenta [3,4]. The conditional DFT wavefunction (2) 
predicts the right amount of vortices inside the quan- 
tum dot but fails to predict the clustering of vortices. 
This can be seen by comparing the ED and DFT re- 
sults for the L — —45 case in Fig. 5. Therefore the L — 
—45 state can be identified as a finite size precursor of 
the v = 1/3 FQHE state only in the ED. 

To conclude, we have studied vortex formation in 
parabolically confined quantum dots and found that 
both the DFT and the ED predict the existence of solu- 
tions with free vortices moving between the electrons. 
We have analysed the structure of these solutions us- 
ing conditional wave functions and we found that the 
DFT gives qualitatively right results but is unable to 
predict the clustering of vortices near the electrons. We 



hope that these results will encourage further analy- 
sis of vortex solutions e.g. in systems of large electron 
number and low symmetry [17]. 
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